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Abstract. An axisymmetric MHD model is examined an- 
alytically to illustrate some key aspects of the physics 
of hot and magnetized outflows which originate in the 
near environment of a central gravitating body. By an- 
alyzing the asymptotical behaviour of the outflows it is 
found that they attain a variety of shapes such as con- 
ical, paraboloidal or cylindrical. However, non cylindri- 
cal asymptotics can be achieved only when the magnetic 
pinching is negligible and the outflow is overpressured 
on its symmetry axis. In cylindrical jet-type asymptotics, 
the outflowing plasma reaches an equilibrium wherein it 
is confined by magnetic forces or gas pressure gradients, 
while it is supported by centrifugal forces or gas pressure 
gradients. In which of the two regimes (with thermal or 
magnetic confinement) a jet can be found depends on the 
efficiency of the central magnetic rotator. The radius and 
terminal speed of the jet are analytically given in terms 
of the variation across the poloidal streamlines of the to- 
tal energy. Large radius of the jet and efficient acceler- 
ation are best obtained when the external confinement 
is provided with comparable contributions by magnetic 
pinching and thermal pressure. In most cases, coUimated 
streamlines undergo oscillations with various wavelengths, 
as also found by other analytical models. Scenarios for the 
evolution of outflows into winds and jets in the different 
confinement regimes are shortly outlined. 

Key words MHD - plasmas - solar wind - stars: mass 
loss - stars: pre-main sequence - ISM: jets and outfiows - 
Galaxies: jets 



1. Introduction 

Nonuniform plasma outfiows seem to be ubiquitous in as- 
trophysics on galactic and extragalactic scales. The clos- 
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est example is the solar wind itself which shows strong 
heliolatitudinal velocity gradients as recently observed by 
Ulysses (Lima & Tsinganos 1996, McComas et al. 1998). 
Further away coUimated outfiows are observed in asso- 
ciation with several galactic objects, such as young and 
evolved stars, planetary nebulae, X-ray binaries and col- 
lapsed objects (for reviews see Ray 1996, Kafatos 1996, 
Mirabel & Rodriguez 1996, Brinkmann & Miiller 1998, 
Livio 1998). Finally, on extragalactic scales jets are ob- 
served to originate in many Active Galactic Nuclei and 
Quasars (Biretta 1996, Ferrari et al. 1996). 

Yet, despite their abundance the basic questions on 
the formation, acceleration and propagation of nonuni- 
form winds and jets have not been fully answered. Never- 
theless, observations seem to indicate that the basic ingre- 
dients for producing astrophysical outflows are some sort 
of heating to launch thermally the wind at the axis plus 
a rotating central gravitating object and/or an accretion 
disk threaded with magnetic fields to accelerate magneto- 
centrifugally and coUimate the outflow. 

1.1. Drivers of the coUimated plasma outflow 

Several mechanisms have been investigated for accelerat- 
ing and coUimating astrophysical outfiows in galactic and 
extragalactic scales. Magnetic rotator forces seem to play 
a rather dominant and crucial role (Lynden-Bell 1996) but 
they are probably not the only relevant mechanism. 
First, thermally driven models are based on the de Laval 
nozzle analogy of the solar wind (Parker 1963, Liffman & 
Siora 1997). This requires the presence of a hot corona 
around the central body of the YSO or the AGN. X-ray 
emission detected in several of these objects may imply 
that thermal effects contribute to the general acceleration 
mechanism at the base of the flow but they are probably 
not the only ingredient. Furthermore, if the wind is asso- 
ciated to a very bright object, the flow can be effectively 
accelerated by the photon flux (radiatively driven winds, 
Cassinelli 1979). Parallely note also that collimation of 
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bipolar outflows from YSOs and Planetary Nebulae by 
external thermal pressure gradients have been extensively 
studied in the frame of the Generalized Wind Blown Bub- 
ble scenario (GWBB, Frank 1998). It has demonstrated 
successfully that magnetic processes may not be the only 
way to achieve collimated outflows. 

Second, magnetic pressure driven models are based on 
the uncoiling spring analogy and have been examined by 
Drains (1983), Uchida & Shibata (1985) and Contopoulos 
(1995). There, it is assumed that a toroidal magnetic field 
i?0 is created and highly amplified by the winding-up of 
its field lines by a radially collapsing and non-Keplerian 
rotating disk. Plasma is then accelerated from the disk in 
the poloidal direction by the action of the resiflting tor- 
sional Alfven waves. 

Third, magnetocentrifugally driven outflow models are 
based on the classical bead on a rotating rigid wire anal- 
ogy. There, the magnetized cold fluid is flung out (even 
to relativistic velocities) from the surface of the Keplerian 
accretion disk, provided that the poloidal field lines are 
inclined enough with respect to the disk axis (Blandford 
& Payne 1982, Pelletier & Pudritz 1992, Contopoulos & 
Lovelace 1994, Cao 1997). This approach is suitable to 
model winds from accretion disks, but is not valid around 
the symmetry axis. Moreover it has been pointed out re- 
cently (Ogilvie & Livio 1998) that, even if the lines are 
sufficiently inclined, a potential barrier still exists that 
can be overcome only by the presence of an extra source 
of energy (e.g. a hot corona). 

In all the above treatments the effects of the combi- 
nation of gas pressure and magnetic fields in accelerating, 
collimating and confining jets have not been discussed ad- 
equately, despite the fact that the role of the gas pressure 
has been recognized for a long time, i.e., that jets are not 
moving in a vacuum (molecular clouds around YSO's, or 
host galaxies in AGN) and hence they must have some 
interaction with the external medium (Ferrari et al. 1996; 
Frank 1998). This approach may also highlight the transi- 
tion from fully thermally driven to fully magnetocentrifu- 
gally driven collimated winds. 

1.2. Radially self-similar models 

As with any fully MHD approach and despite of the simpli- 
fications of steadiness and axisymmetric geometry, several 
approximations are still unavoidable in order to obtain 
exact solutions useful for an understanding of the MHD 
mechanism for the initial acceleration and final coUima- 
tion. Thus, one simple analytical way out is the use of 
self-similarity This hypothesis allows an analysis in a 2- 
D geometry of the MHD equations which reduce then to 
a system of ordinary differential equations. The basis of 
the self-similarity treatment is the assumption of a scaling 
law of one of the variables as function of one of the coor- 
dinates. The choice of the scaling variable depends on the 
specific astrophysical problem. 



Several models self-similar in the radial direction have 
been investigated to analyze the structure of winds from 
accretion disks (Blandford & Payne 1982, Contopoulos & 
Lovelace 1994, Li et al. 1992, Li 1995, 1996, Ferreira 1997, 
Ostriker 1997). In these models the driving force and the 
collimation derive from a combination of the magnetic 
and centrifugal forces. Moreover, as disc-winds arc associ- 
ated with jets, these studies usually do not consider under 
which parametric conditions full collimation is obtained. 
Exceptions are given in Pelletier & Pudritz (1992) and 
Contopoulos & Lovelace (1994) where the collimation effi- 
ciency is linked to a current fiowing in agreement with the 
Heyvaerts & Norman (1989) general analysis. However, 
the absence of an exact crossing of all the existing critical 
points in the solutions presented in these papers prevents 
from considering their conclusions as definitive. Neverthe- 
less, it has been shown that within the frame of self-similar 
disc-wind assumptions, it is possible to cross all critical 
points thus getting meaningful solutions (Tsinganos et al. 
1996, Vlahakis 1998). Moreover, the role of the inhomo- 
geneity in the pressure distribution has not been taken 
into account until recently in these models (Ferreira 1997) 
and a full parametric study of this extra variable is yet to 
be performed. 

l.S. Meridionally self-similar models 

In a series of studies, solutions of the MHD equations that 
are self-similar in the meridional direction have been also 
analyzed (Tsinganos & Trussoni 1990, 1991, Tsinganos & 
Sauty 1992a, b. Papers I and II of this series, Trussoni 
& Tsinganos 1993, Sauty & Tsinganos 1994, Paper HI 
of this series, Trussoni et al. 1997, henceforth TTS97). 
Such a treatment allows to study the physical properties 
of the outflow close to its rotational axis. As in this region 
the contribution to acceleration of the magnetocentrifu- 
gal forces is small, the effect of a thermal driving force is 
essential. This implies also that the structure of the gas 
pressure in the flow is essential. 

Two main classes of such self similar solutions have 
been investigated depending on whether the components 
of the pressure gradient along the radial and meridional 
directions are related or not. In the second case the shape 
of the streamlines and fieldlines is prescribed 'a priori', 
and the main features of the dynamical variables are self- 
consistently deduced from the integration. In particular, it 
has been shown that acceptable solutions for magnetized 
flows with asymptotic superAlfvenic velocity exist only 
when rotation is included (Tsinganos & Trussoni 1991, 
Trussoni & Tsinganos 1993, TTS97). As a consequence of 
this study it seems that even pressure confined jets from 
slow magnetic rotators need magnetic fields and rotation. 

In the other case, in which the two components of the 
gas pressure are related, the structure of the streamlines 
is deduced as a self-consistent solution of the MHD equa- 
tions. It has been shown (Papers I and II) that hydro- 
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dynamical and nonrotating magnetized winds are always 
radially expanding from the source. On the other hand, 
rotating magnetized flows with a spherically symmetric 
structure for the pressure gradient can have final super- 
Alfvenic velocities with either radial or coUimated asymp- 
totically streamlines, depending on the values of the pa- 
rameters (Paper III). This allows to deduce a criterium to 
select conically expanding winds from cylindrically coUi- 
mated jets (Sauty et al. 1996). 

1.4- Plan of this paper 

We extend here the analysis of Paper III to the more gen- 
eral case of solutions for rotating magnetized winds with a 
nonspherically symmetric gas pressure. In the present pa- 
per we concentrate on the asymptotic analysis and its link 
to the initial boundary conditions: this allows us to derive 
a general criterion for the coUimation of winds into jets. 
The analysis of the properties of the complete numerical 
solutions deserves a separate study which is postponed to 
a following paper. 

In Sec. 2 we summarize the properties of the merid- 
ional self-similar MHD equations while in Sec. 3 we dis- 
cuss the energetic structure of the outflow. In particular we 
show that an energy integral exists that links the asymp- 
totic regime to the boundary conditions at the base, allow- 
ing to formulate a general criterion for the collimation of 
the wind. The different physical conditions for asymptotic 
confinement (magnetic or thermal) are discussed in detail 
in Sec. 4, and in Sec. 5 we show that oscillating configu- 
rations can be present in cylindrically coUimated jets. In 
Sec. 6 the equilibrium asymptotic properties of non coUi- 
mated flows are outlined, while in Sec. 7 we summarize the 
results and shortly discuss the astrophysical implications 
of our analysis. 

2. Meridionally self-similar MHD model 

2.1. Steady axisymmetric ideal MHD outflows 

The global dynamical properties of cosmic winds and jets 
are usually analyzed by assuming that they represent out- 
Uows of a fully ionized plasma with a bulk speed V and 
carrying a magnetic field B in the gravitational field of a 
central body of mass Ai . The familiar MHD equations are 
employed for a physical description of these phenomena. 
In particular, under steady and axisymmetric conditions 
{d/dt = d/dtp = 0), the MHD equations are known to ad- 
mit certain free integrals, i.e., functions which remain con- 
stant on the magnetic surfaces generated by the revolution 
around the magnetic/flow symmetry axis of the system of 
a poloidal magnetic line A{r^ 9) = constant (Tsinganos 
1982). Specifically, on the surface of such a flux tube A = 
const., the following physical quantities remain invariant 
throughout the extent of these surfaces from the base to 
infinity: 

— ^a(A), the ratio of the magnetic and mass fluxes. 



— L{A), the total specific angular momentum carried by 
the flow and the magnetic field, 

— ^(A), the corotation frequency or angular velocity of 
each streamline at the base of the flow. 

Furthermore, it is well known that the poloidal {p) 
and azimuthal (toroidal, ip) components of the magnetic 
field and the velocity can be expressed in terms of these 
free integrals and the poloidal Alfven Mach number, using 
spherical (r, 9, ip) or cylindrical {w, ip, z) coordinates (for 
details see Paper III). In particular, the poloidal Alfven 
Mach number (or Alfven number) M is. 



K' 



^' 



' Bl Airp 



(2.1) 



On the other hand, the two integrals L{A) and il(A) are 
not independent ij the flow is transalfvenic. In such a case, 
at the cylindrical distance vja of the Alfven point (M = 1) 
from the field/flow axis of a flux tube labeled by A they 
are related as vj1^{A) = (r, sin^a)^ = L/U. 

2.2. Generalized Bernoulli integral 

A fourth constant of the motion expresses the conservation 
of energy along streamlines. Thus, by projecting the mo- 
mentum equation along a streamline, taking into account 
the first law of thermodynamics for energy conservation, 
we obtain the generalized classical Bernoulli integral (Pa- 
per HI), 



2 ■ P 2 f 



where 

e;„(A) 



—ruB^ + h-e:^{A), 
(2.2a) 

(2.2b) 



qir',A) 



dr' 

,^p{r',A)Vr{r',A) ' 

h is the enthalpy of the perfect monoatomic gas (F — 5/3), 
q is the net local volumetric heating/cooling rate, and Q 
the gravitational constant. Thus, at a given radial dis- 
tance r along the streamline labeled by A, the conserved 
energy E{A) represents the sum of the kinetic, gravita- 
tional, Poynting and ideal thermal energy flux densities 
per unit of mass flux density, minus the extra heat re- 
ceived by the flow between the anchored footpoint at a 
basal radial distance Tq and the point r under considera- 
tion, e;^(A). 

2.3. Self- similarity: scaling laws for the variables 

The model analysed in this paper belongs to the wide 
class of meridionally self-similar MHD equilibria (see also 
Trussoni et al. 1996; Tsinganos et al. 1996; TTS97; Vla- 
hakis & Tsinganos 1998, henceforth VT98). In the follow- 
ing we briefly summarize the main steps for the construc- 
tion of such a model (see Appendix A for more technical 
details). 
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For convenience, first of all the variables are normal- 
ized to their respective values at the Alfvcn surface along 
the axis of rotation, r = r* . In particular, we define the di- 
mensionless radial distance R ~ r/r^, and the Alfven speed 
y* = Bl/Anp^, where S*, K and p* are the poloidal mag- 
netic field, poloidal velocity and density along the polar 
axis at the characteristic radius r^ . For the magnetic flux 
function A we define its dimensionless form by 



a{R,( 



A{r,e) 



Note that along the polar axis a{R, 0) — 0. To obtain the 
final expressions for the physical variables, we make the 
following crucial assumptions: 

— First, we assume that the Alfven surface is spherical, 
M = M{R). Then, according to Eq. (2.1), the density 
can be expressed as the product of a function of R [i.e. 
1/AttM^{R)] and a function of a [i.e. '^\{a)]. Further- 
more, we Taylor expand the function ^^(a) to the first 
order in a such that the variation of the density on a 
spherical surface of given radius R is proportional to 
the magnetic flux a. 

— Second, we assume that the magnetic flux function a is 
expressed as the product of a function of R and a func- 
tion of 0. Furthermore, for the function of 9 we take 
a dipolar dependence with the colatitude 6. This im- 
mediately implies that the Alfven cross sectional area 
TTZu'^ of a flux surface is proportional to the correspond- 
ing magnetic flux a. Also, the ratio G^ = vd^ lvo\ of the 
cross sectional area of the flux tube to the Alfven cross 
sectional area of the same flux tube depends solely on 
the radial distance R. 

— Third, we assume that the total axial current Iz en- 
closed by a flux tube a = const, is proportional to the 
corresponding magnetic flux. This assumption fixes the 
angular momentum integral L (Paper III). Note that 
at once the integral of the corotation frequency 51 fol- 
lows from its relation with L at the Alfven distance, 
L = flru'^ . Note also that the integrals L and il are cho- 
sen such that L'^A and L0\['^ contain only first order 
a-terms, in analogy with the previous assumptions. 

— Fourth, we assume that the a-dependence of the gas 
pressure is similar to that of the density distribution. 
This means that the pressure is ultimately a function 
of the density along a given magnetic surface, a situ- 
ation analogous to the often used polytropic assump- 
tion. However, this implicit relationship between pres- 
sure and density is much more general than the some- 
how artificial polytropic assumption. Contrary to the 
polytropic relation, its exact form is not imposed a 
priori but is determined by the full solution. 

Altogether, the four main assumptions of this merid- 
ionally self-similar model can be summarized as follows. 



p{R,a) 



P* 



NP{R) 



(l + (5a) 



*i 



47rp,(l + (5a), (2.3a) 



vj'^{R,a)^rlG^{R)a, zul{a) ^ r^a , (2.3b) 

L*A = Ar,B,a , Ln^\ = X^B^a , (2.3c) 

PiR, a) = ip,K'n(i?)(l + Ka) . (2.3d) 

The introduced parameters 6, n and A measure the vari- 
ation with the colatitude of the density, pressure and ro- 
tation, respectively. A fourth parameter v enters from the 
momentum equation as the ratio, at the Alfven distance 
along the polar axis, of the escape speed to the flow speed 
there, 

2 _ 2gM 



r.Vi 



(2.4) 



2. 4- Magnetic rotator energy 



An important physical quantity in magnetized outflows 
is the so called magnetic rotator energy (Michel 1969, 
Belcher & McGregor 1976), 



E 



MR 



= 9.L. 



(2.5a) 



The basal Poynting energy i?poynt..o, deflned as the ratio of 
the Poynting flux density S^ per unit of mass flux density 
pVz, is roughly equal to the magnetic rotator energy flL 
if at the base the radius of the jet is much smaller than 
the Alfven radius (Go ^1) and the Alfven number there 
is also negligibly small [Mo <C 1), 



-L/Poynt.,o — T^ 



Sz 
pVz 



ni 



1 - M 2 



VLL. 



(2.5b) 



Let Eo be the sum of the kinetic, gravitational and thermal 
energies per unit mass at the base of the outflow. Then 
the total available energy for the outflow at the base is 
E K, Eo + ^L. Accordingly, we have an outflow from a Fast 
Magnetic Rotator (FMR) when Eo ^ JIL and an outflow 
from a Slow Magnetic Rotator (SMR) in the opposite case 
of Eo > riL. 

2.5. Solving the self-similar MHD equations 

In order to solve the resulting MHD equations, it is useful 
to introduce an extra function F(R) (Papers II and HI), 



F{R) = 21 



din 



G{R) \ 



d Ini? / 



(2.6) 



Evidently, while G{R) defined in Eq. (2.3b) measures 
the dimensionless cylindrical radius of a flux tube at the 
distance R, F{R) is simply giving the expansion factor 
of the streamlines. The limiting case F{R) — corre- 
sponds to conical expansion and radial fleldlines, while 
for F{R) = 2 we have cylindrical expansion parallel to the 
axis (collimation). In between these two regimes the flow 
is paraboloidal. 

The above assumptions, Eqs. (2.3), immediately give 
the components of the velocity and magnetic flelds (Eqs. 
A. 3 in Appendix A). On the other hand, the momentum 
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conservation law in combination with the above assump- 
tions gives four ordinary differential equations for the four 
variables M'^{R), F{R), n(i?) and G{R) (see Appendix A 
for details). 

The complete solution of these equations, from the 
base of the outflow to infinity, with the required crossing of 
all appropriate critical points, is indeed an interesting un- 
dertaking and worth of a separate paper. Here instead, we 
shall concentrate on some novel results obtained solely by 
solving the equations asymptotically far from the Alfvcn 
surface (R ^ 1) and taking into account the boundary 
conditions on the source. 



It is worth to digress for a moment and try to get some 
insight into the physical meaning of these two conserved 
components of the specific energy, £ and AS. 

3.1.1. Polar specific energy 

In the first expression, Eq. (3.3a), the polar energy flux £ is 
composed of four terms which are successively the poloidal 
(i.e. radial here) kinetic and gravitational energies, the 
enthalpy and the heating along the polar axis. 

The polar specific energy £ can be evaluated at both 
the base of the wind Rq and far from it as R — > oo, 



3. The energy integral and collimation criterion 

3.1. The generalized Bernoulli integral 

A nonadiabatic fiow of a monoatomic gas with ratio of 
specific heats F — 5/3 is always heated at a net volumetric 
rate q 



q^ pV ■ [\/h 



VP 

P 



r p 

f^7 



(3.1a) 



With expressions (2.3a) and (2.3d) for the gas density and 
pressure, it follows immediately that this heating can be 
written as (see Sec. 5 in Paper III for details of the deriva- 
tion), 



q{R,a) 



KM 



p{R,a)Vr{R,a) 



2r* 1 + Sa 



Q{R) 



(3.1b) 



where the dimensionless specific heating rate per unit of 
radial length along a given streamline is. 
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' dn 

^^ -TT^ + rn- 
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dR 



(3.1c) 



Hence, the generalized classical Bernoulli integral (2.2) 
takes the simpler form 



^ ' 2 * l + Sa 



(3.2) 



where the two constants £ and A£ represent the polar 
specific energy and the variation across a streamline of the 
specific energy, respectively (in Paper HI £ was denoted 
by .fi, AS by .^2 and Q by Qi). 

It is straightforward to show from Eqs. (2.2), (2.3), 
(A. 3) and (3.1)-(3.2) that £ and A£ have the following 
analytical expressions (Paper III) 



£ = 
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Q{R)dR, (3.3a) 
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Ro 



T 



T - 1 



^oMl 



(3.4a) 
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Q{R)dR. (3.4b) 



Ro 



At the base, wherein the kinetic energy of the outflow is 
negligible, Eq. (3.4a) shows that the polar energy has ba- 
sically two terms: the gravitational energy and the initial 
input of thermal energy in the form of enthalpy. On the 
other hand, at infinity, Eq. (3.4b), the conserved polar spe- 
cific energy is composed of the final kinetic energy along 
the polar axis and the terminal enthalpy minus the addi- 
tional extra heating which the flow has received during its 
propagation from Ro to infinity. 

Note that if the wind is cylindrically coUimated, Moo, 
Goo and Hoo have finite values. In all other cases. Moo and 
Goo rnay be unbounded, although their ratio, which is the 
polar speed in units of the Alfven speed, should remain 
finite, 

^- ^- (3.4c) 



r^2 



K 



Moreover the terminal pressure Hoo vanishes unless the 
integral of the heating diverges, a rather unphysical situ- 
ation corresponding to an infinite input of heat. 

The conservation of the polar energy simply expresses 
the fact that the flow along the polar axis is thermally 
driven. Furthermore, from Eqs. (3.4a,b) it becomes evi- 
dent how the conversion of the heat content of the plasma 
into kinetic and gravitational energy maintains the out- 
flow, 



-(n^M^-HooM^ 



Q{R)dR 



v^ M4 

\ — 

^o tj^oo 



(3.5) 



In other words, the decrease of the enthalpy at infinity to- 
gether with the external heat input integrated along the 
polar streamline, on one hand lifts the gas out of the grav- 
itational potential well and on the other, gives to it a finite 
terminal speed. Of course, this is nothing more than the 
classical picture of the Parker thermally driven wind. 



Sauty et al. Nonradial and nonpolytropic astrophysical outflows. IV. 



3.1.2. Variation of the specific energy across streamlines 

Tlie second conserved component A£ of the specific en- 
ergy gives the excess or deficit of the volumetric total en- 
ergy E at a nonpolar streamline as compared to the corre- 
sponding energy at the polar axis and the same spherical 
distance, normalized to the volumetric energy of the mag- 
netic rotator. Thus, AS has five contributions which cor- 
respond to the five different terms appearing successively 
in the RHS of Eq. (3.3b). Each one represents the varia- 
tion - in units of the volumetric energy of the magnetic 
rotator - between any streamline and the polar axis of 
(i) the poloidal kinetic energy, (ii) the volumetric gravita- 
tional energy, (iii) the azimuthal kinetic energy (which is 
zero along the polar axis), (iv) the Poynting flux (which is 
also zero along the polar axis) and (v) the thermal content 
(enthalpy plus heating; see Appendix B for details). 
In a more compact way we may write AS as follows. 



AS _ p{R,a)E{a 
2A2 " 



/9(i?, pole)i?(pole) 



p{R, a)L{a)fl(a) 



(3.6) 



Evidently, AS/2X'^ represents the variation across the flow 
of the total volumetric energy in units of the volumetric 
energy of the magnetic rotator. Therefore, the sign of AS 
determines whether there is a deficit of energy per unit 
volume (and not per unit mass) along the polar streamline 
as compared to the other streamlines (case AS > 0) or an 
excess of energy in the polar streamline as compared to 
the other nonpolar streamlines (case AS < 0). 

Furthermore, AS /2X^ can be expressed (see Appendix 
B for more comments) in terms of the conditions at the 
source boundary Ro where the cylindrical radius is Wo {a) , 
the escape speed Vesc,o, the polar density po(pole) and the 
density at any other streamline po{a)'- 

AS _A [po{Epoynt.,o + En.o + ho + Eg,o)] .„ „ ^ 

— 2 ~ aTTI? ^ ' (,o./aj 



2A 



A{pEmr)o 



where A{pEmr)o is the variation of the energy of the mag- 
netic rotator, A(/?i?Poynt.)o is the variation of the Poynting 
energy, A{pEYi)o is the variation of the rotational energy 
at the base, A(p_Eg)o is the variation of the volumetric 
gravitational energy at the base and A{ph)o is the vari- 
ation of the volumetric thermal flux at the base, respec- 
tively, 

A(p£;poynt)o = Po{a)Epoynt.oi<^) = Po{a)i^L - il'^zul) , 

(3.7b) 

A{pEji)o = Po(a)£^R,o(a) = Po(a)^^ = Po{a)—^ — , 



A{pEq)o = [Po{a) - po(pole) 



(3.7c) 
(3.7d) 



Mph)o = p^[^o(a) - Poipolc)] , (3.7e) 

A{pEmr)o = po{a)EMR{a) = po{a)L{a)Q{a) . (3.7f) 



In this notation, A always denotes a variation across the 
fieldlines at a given radial distance R, i.e. Ay = y{R, a) — 
y{R, pole) for every function y{R, a). 

In Eq. (3.7a) note that (see also Eqs. 2.5) 



-Epoynt.,o + Ep^o — Emr — Ep,o > . 



(3.8) 



The Poynting flux plus the rotational energy is simply 
the energy of the magnetic rotator minus the rotational 
energy. This last form is the one used in Paper III. In other 
words, and even in the slow magnetic rotator limit, the 
rotational energy never exceeds the energy of the magnetic 
rotator. 

3.2. Energetic definition of Efficient/Inefficient Magnetic 
Rotators 

At this point we inevitably note that S and AS are two 
inconvenient constants because their absolute values de- 
pend on the integration of the total heating supply and 
so they can be evaluated only after the problem has been 
solved and the required heating can be calculated. How- 
ever these two constants are related to each other. In fact, 
the last two terms in the expression of AS in Eq. (3.3b), 
which correspond to the transverse variations of enthalpy 
and heating, are identical to the last two terms of S within 
a factor of k. Evidently, this is due to the assumptions on 
the pressure and density distribution, Eqs. (2.3a,d). These 
initial assumptions imply the existence of an implicit re- 
lationship between the latitudinally normalized pressure 
and density, 



P{R,a) 



1 



p{R, a) 

p{R,o) 



1 



(3.9) 



The situation is akin to the more familiar polytropic 
ansatz, although there the relationship between pressure 
and density is explicit. In TTS97 the generalized Bernoulli 
integral has indeed a form similar to Eq. (3.2), but the two 
constants S and AS are not related to each other as in Eqs. 
(3.3), because the spherically symmetric part of the pres- 
sure is not related to the corresponding nonspherical part. 
For this reason, it was impossible to find a relationship be- 
tween P and p of the form of Eq. (3.9) and therefore any 
convenient form of the Bernoulli integral. 

With this in mind, we can eliminate from the expres- 
sions of S and AS in (Eqs. 3.3) the inconvenient enthalpy 
and heating terms (Paper III) by defining the new con- 
stant 

e = AS-KS. (3.10) 

Now this quantity e, in addition of being a constant for 
all streamlines. 



M" 



[GRy 



4 



A2 fM^-G^ 



G2 V 1 - M2 



2X' 



'~G^ 
1 



{5 - K)v'^ 



R 



G' 



1 -M2 



(3.11) 
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can be calculated a priori from the conditions at the base 
of the outflow, without a need to know the total input of 
heating along each line. 

A careful look at Eq. (3.11) shows that all the trans- 
verse variations of the total energy, simply reproducing, 
within a scaling factor k, the effect of thermally driven 
winds along the pole (Eq. 3.4), have been removed (see 
Eq. B.7 in Appendix B). 

In fact, comparing Eq. (3.11) to Eq. (3.3b), we see 
that e contains the same terms as A£ except the heat 
content, but with two extra terms proportional to n. The 
first of these two terms {kM^ /G^) represents simply the 
transverse variation of the heat content which is converted 
into kinetic energy in a thermally driven wind, as seen by 
Eq. (3.5). The second term (kv"^ / R) is the variation with 
the latitude of the thermal energy which along the pole 
supports the plasma against gravity. 

Let's assume for a moment (5 = k > 0, such that the en- 
thalpy and the temperature (ex P/ p) are spherically sym- 
metric. Since the pressure is larger on a nonpolar stream- 
line, we have higher heating rate q there: the extra heating 
converted into kinetic energy is kM^ jG^ (Eq. 3.5). In the 
total energy variation budget it represents the efficiency of 
thermal confinement. Therefore it must be removed from 
the energy variation in order to form the constant e. The 
same holds if k < except that this k term will tend to 
decollimate the outflow. 

Now, if K = and (5 > we see that there is an excess 
of gravitational potential —5v^/R because the plasma is 
heavier on a nonpolar streamline. In order to achieve equi- 
librium, part of the Poynting flux and part of the centrifu- 
gal energy must compensate this term. This reduces the 
energy available for magnetic confinement. If i5 > k, we 
need to correct the previous argument because part of the 
weight of the plasma is supported on a non polar stream- 
line by an increase of the pressure gradient. This compen- 
sation is exactly kv^ j R. Thus the term —{5—k)v'^/R is the 
effective increase of the gravitational potential that must 
be compensated by some non thermal drivers, the mag- 
netic driver for instance. It reduces the efficiency of the 
magnetic rotator to coUimate the flow. Similar arguments 
hold if (5 < or (5 < K. 

As in Eqs. (3.7) let's express e/2A^ in terms of the 
conditions at the source boundary Tq (assuming again that 
the velocity is negligible there, see Appendix B for details 
of the derivation) , 

where iSpoynt. and E^o have been already defined. A£^q 
is the excess or the deficit on a nonpolar streamline com- 
pared to the polar one of the gravitational energy (per unit 
mass) which is not compensated by the thermal driving. 



l^E*r 



GM 



Tpja) 
To (pole) 



GM (S - K)a 



Sa 
(3.12b) 



It is indeed the term proportional to {S — k)u'^ in Eq. (3.11) 
and the symbol A keeps the same meaning as previously 
(see Appendix B). 

It is worth to remark that this corrected gravitational 
term plays an important role in thermally accelerating the 
flow (Tsinganos & Vlastou 1988; Paper I) because it is pro- 
portional to the relative variation of the temperature. We 
know from previous numerical studies that (k — S) ought 
to be negative in order that we have efficient initial accel- 
eration along the polar axis. This amounts to say that the 
temperature along the polar axis must be larger than the 
temperature along a non polar line. Then, the corrected 
gravitational term in Eq. (3.12a) is always negative such 
that it must be compensated with part of the initial input 
of the magnetocentrifugal terms (Poynting and rotational) 
at the base of the flow. 

Hence, e > means that the magnetocentrifugal terms 
are dominating the variation of gravity and that there is 
some energy left from the magnetic rotator to coUimate 
the wind. While e < means that the magnetic rotator 
cannot coUimate the wind by itself. Of course the coUimat- 
ing efficiency of the magnetic rotator may be eventually 
lowered if there is further pressure gradient acting out- 
wards in the wind [k < 0) but e/2A^ really quantifies the 
original strength of the magnetic rotator to support the 
coUimation of the fiow. 

As a conclusion of this subsection, we may define as 
Inefficient Magnetic Rotators (IMR) the magnetic 
rotators which are not able to confine the fiow through 
magnetic processes alone and have e < 0. Conversely we 
shaU caU Efficient Magnetic Rotators (EMR) the 
magnetic rotators potentially able to magnetically confine 
the fiow and which have e > 0. We shall further illustrate 
this definition at the end of the next subsection. Within 
this definition the classical Slow Magnetic Rotators 
(SMR) and Fast Magnetic Rotators (FMR) corre- 
spond respectively to (IMR) and (EMR) but only in the 
limit where all other energies are distributed in a spheri- 
cally symmetric manner at the source base. 

3.3. Energetic criterion for cylindrical coUimation 

The coUimation of an outflow can be either of magnetic, 
or of thermal origin. In the following, we discuss how to 
measure the distribution of the thermal content along and 
across the flow, before reaching some conclusions on the 
coUimation itself. 

In a thermally driven wind, all thermal input (internal 
enthalpy plus external heating provided along the flow) is 
not necessarily fully converted into other forms of energy, 
unless the terminal temperature is exactly zero. There 
always remains some asymptotic thermal content in the 
form of enthalpy. Conversely, we can deflne the heat con- 
tent that is really used by the flow by defining the con- 
verted enthalpy 



Mr, A) = h{r, A) + S^iA) - h{oo, A) 



(3.13a) 
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Along a fieldline A and sonic radius r, the converted en- 
thalpy h is simply the sum of the enthalpy of the gas at 
this point and the external heat which it will receive on 
its way to infinity, Q'^{A) (Eq. 2.2b), minus the enthalpy 
that will still remain in the gas asymptotically. Note that 
in the polytropic case this converted enthalpy is simply the 
variation along the flow of the effective enthalpy ho — hoo ■, 
where the adiabatic index F is replaced by some effective 
7 < r, as explained in Paper III. Then we can define a 
constant along each streamline 

E{A) = \v^ + \v'^-^-^-^^r.B, + h. (3.13b) 

We may also redefine the variation across fieldlines of the 
volumetric energy normalized with the energy of the mag- 
netic rotator, but including the converted enthalpy which 
will indeed be used by the flow, instead of the enthalpy. 
In other words, wc may define a new quantity e' in the 
way we defined A£ in Eq. (3.7a), but using the converted 
enthalpy instead of the enthalpy. 



e' _ p{R, a)E{a) - p{R, polc)i?(polc) 
2A2 ~ p{R, a)L{a)n{a) 

Thus we have at the base 

g/ A Po{Epoynt.,o + Ej\,o + ho + Eq^o) 

2A2 " 



A{pEMVi)o 



(3.14a) 



(3.14b) 



where all the terms have the same meaning as in Eq. (3.7a) 
except the transverse variation of the total converted en- 
thalpy of the flow which is simply 

A{ph)o == po{a)ho{a) - po (pole) /iq (pole) . (3.14c) 

Working out this definition together with Eqs. (3.7), (3.3) 
and (B.5.), wc find the following relation 



e' == A£ - 



r 



-n^M^ + K Q{R)dR. (3.15) 
^ Jro 



Thus e' is simply the difference of A£ and the total heat 
content of the flow at infinity. As a consequence, in our 
model e' is a constant which can be evaluated a priori 
at any r using Eq. (3.15). In the particular case of a flow 
which is asymptotically cylindrically coUimated {Foo = 2), 
this parameter can be evaluated at infinity in a simple way 
(cf. Eqs. 3.3b -3.15), 



{Ml 



G\ 



2 \2 



Gl 



1 



2A2 2GL(M^-1)2 • Ml-l 



(3.16) 



Note that now in e' there are left only the magnetocen- 
trifugal terms: variation of the azimuthal kinetic energy 
and Poynting flux. The important result is that this new 
parameter is always positive in a cylindrical jet, if the jet is 
asymptotically super Alfvenic, Moo > 1, and transalfvenic 
so the asymptotic radius is larger than the Alfvcn radius. 



Goo > 1. In other words, a necessary condition for achiev- 
ing cyhndrical collimation is that e' > 0. 

The criterion for cylindrical collimation is thus explic- 
itly equivalent to the criterion given in Paper III, ex- 
cept that now we have included the thermal contributions: 
cylindrical collimation can be achieved only if there is an 
excess of energy on a non polar line compared to the polar 
one. However, it is not the variation across the lines of the 
total thermal energy input that enters in the definition 
of the criterion, but the variation of the thermal energy 
that is effectively converted into some other form of energy 
between the base and the asymptotics (Eq. 3.13). 

Two contributions may arise to give a positive value 
for e': either because the energy of the magnetic rotator 
dominates as in Paper III, or because the thermal con- 
tribution converted to non thermal energy in the flow is 
higher outside the polar axis. This last point may be bet- 
ter realized if we note that 



* oo 



(3.17) 



Thus, e' splits into two parts. The first is e which is es- 
sentially positive when the energy of the magnetic rotator 
dominates (see Paper III and the previous subsection). 
The second corresponds to the variation with colatitude 
of the thermal energy that has been converted to kinetic 
energy (see Eqs. 3.5 and 3.4c). 

Altogether, there are two ways to have e' > 0: 
Either, when the outflow is magnetically dominated, 
which means that e is positive and k may be either positive 
(which adds some extra pressure confinement) , or negative 
(which corresponds to pressure support of the jet) within 
some limits. 

Or, conversely, when there is a significant contribution of 
the variation of the enthalpy-|-heating term that is con- 
verted into kinetic poloidal energy, then kV^/V^ is large 
enough which implies k > while e may be negative. This 
does not necessarily implies that the flow is pressure con- 
fined as we shall see later. 

4. Asymptotic confinement 

We proceed now towards an asymptotic analysis of the 
equations of motion, in particular in the case of cylindrical 
collimation. 

4-.1. Asymptotic equilibrium in cylindrically coUimated 
outflows 

When the equilibrium is asymptotically cylindrical e' > 
as discussed above. Taking the dominant terms in the 0- 
component of the momentum equation or, equivalently by 
expressing force balance across the cylindrical fieldlines, 
we obtain the condition of MHD equilibrium in the cylin- 
drical radius direction w expressed by the equation. 



fB + fp = 0. 



(4.1) 



c)n =0.1 

•* on 



10" 



0.1 
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Fig. 1. Plots of the asymptotic cylindrical radius normalized to the cylindrical radius at the Alfven surface Goo vs. e/2A for 
various values of the flnal pressure: IIoo = (panel a), = 0.01 (panel b), = 0.1 (panel c) and = 1 (panel d). Each curve is 
drawn for a constant value of k/2A^ between —0.1 and 0.1 which labels the curve. On the left of the dotted line is the domain 
of pressure confined jets while on the right of the dashed line is the domain of pressure supported jets and in between is the 
domain of magnetocentrifugal jets. 



In the asymptotic regime, the centrifugal (/c)j magnetic 
{f b) and gas pressure gradient [f p) forces have the fa- 
miliar expressions, 



(4.2a) 
(4.2b) 
(4.2c) 



fc 


= 




d / 

dm \ 


Sit 


\ R2 1 

/ Attw 


fp 


= - 


dP . 

dw 



w . 



In our notation, they can be written as, 

r2 \2 / 7i#2 ,-i2 \ 2 



Sc = 



p.vi y? (Ml 



Gl 



Gl,Ml \ Ml-1 



p.V^2X^ fCl-l 



G3 I ]\,p _ I 



(4.3a) 



(4.3b) 






/azu . 



(4.3c) 



Note that always the centrifugal force acts outwards while 
the total magnetic force (pinching plus pressure) inwards. 
On the other hand, the last term (not appearing in the 
Paper III study) is the pressure gradient that acts out- 
wards if the flow is overpressured (k < 0, i.e. the pres- 
sure decreases away from the axis). In this case the jet is 
necessarily magnetically confined but either centrifugally 
supported or pressure supported. Conversely, the pressure 
gradient acts inwards if the flow is underpressured (k > 0, 
i.e. the pressure increases away from the axis). In this case 
the flow is centrifugally supported but may be either mag- 
netically confined or pressure confined. 
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Fig. 2. Plots of the asymptotic Alfven number Moo vs. e/2A^ for the same values of IIoo and k/2A^ as in Fig. 1. On the left of 
the dotted line is the domain of pressure confined jets while on the right of the dashed line is the domain of pressure supported 
jets, in between is the domain of magnetocentrifugal jets (except for a small overlap around the dashed line of panel b) . 



By combining the asymptotic transverse force balance 
(Eqs. 4.1 - 4.3), 



1 



2G^M^ 



M'i 



Gl 



Ml, 






1 



2A2 



n„o = o, 



(4.4a) 



with the expression of e calculated at infinity (Eq. 3.11) 



e 
2A2 



n Mt 



1 



2x^Gt [Mi^iy 



{Ml-Gl^f 



2Lr^, 



+ {Gl,~l)(Ml~l) 



(4.4b) 



we obtain the asymptotic jet radius and Alfven number as 
functions of the parameters e/2X^, k/2X^ and the asymp- 
totic pressure IIoo. Plots of the resulting values of Goo, 
Moo and the axial terminal speed T^oo/K vs. e/2\^ for 
four representative values of IIoo (0, 0.01, 0.1 and 1) are 
shown in Figs. 1-3. In each of these panels the values of 



Goo , Moo and Foo are plotted for a range of values of the 



pressure parameter k/2A2 between 
label the curves. 



-0.1 and 0.1 which 



In parallel we also plot (Fig. 4) Goo, Moo and Foo/K 
vs. e' /2X'^ using Eq. (3.17) to determine e' /2X^. Note that 
(see Eqs. 4.4) Goo, Moo and Ko/K depend only on e' /2X^ 
and Knoo/2A2. So conversely to Figs. 1-3 where each curve 
is drawn for given values of IIoo and k,I2X^ independently, 
the curves of Fig. 4 are drawn for a constant and unique 
value of kIIoo/SA^ which labels the curve. In Fig. 5 we 
make an explicit comparison of the plot Goo vs. e/2X'^ 
of Fig. la and the corresponding curve Goo{^' /2X^) for 
k\Ioo/2X^ = 0. 

In these plots we may find three different regimes for 
the asymptotic state of the collimated outflow according 
to the various confinement and support conditions across 
the jet. 
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Fig. 3. Plots of the asymptotic axial velocity normalized to the polar Alfven velocity Voo/V, vs. e/2\^ for the same values of IIoo 
and k/2A^ as in Fig. 1. The domain of pressure conflned jets is clearly on the left side of the dotted line and the magnetocentrifugal 
one in between the dotted and dashed line. However the pressure supported domain overlaps the magnetocentrifugal one. 



4-. 2. Magnetocentrifugal jets with fc = |/b| 

This is the case when the pressure gradient is exactly zero 
fp = i.e., the pinching magnetic force is balanced by the 
inertial (centrifugal) force alone. The two following cases 
correspond to this situation. 



Emax the jet's radius, Alfven number and terminal speed 
increase. Finally, as e — > 0, Goo — > oo, Mao — * oo 
and Voo — > cx). The streamlines become conical and the 
asymptotic speed diverges. No collimated solutions at all 
exist for e < 0, as is evident from the plots. 



4.2.1. Spherically symmetric pressure, k — 

The pressure is everywhere spherically symmetric. The 
terminal value of the pressure does not affect the asymp- 
totic equilibrium in the jet, regardless if it is finite (Hoc J^ 
0) or zero (IIoo = 0). This situation corresponds to the 
thick solid curve labeled in Figs. 1-3. This special 
and simplest case has been already discussed in detail 
in Paper III where it was found that the outflow coUi- 
mates into a cylindrical jet only for e > 0, because in 
this case e = e' 
perAlfvenic, M. 
e/2A2 — > emax/2A2 = (2 - ^/2) « 0.586, then Afoo — > 1 
and G^ — > 1 (see Figs. 1 and 2). As e decreases from 



' For a given A, if the flow remains su- 
> 1, an upper limit exists for e. For 



4.2.2. Vanishing asymptotic pressure, IIoo = 

The asymptotic gas pressure is zero such that we have 
again fp — 0, even though k y^ 0. This is the case shown 
in Figs. la-3a where, besides the collimated solutions of 
the Paper III case obtained for e > and k = (the thick 
solid branch), we have now collimated solutions for prac- 
tically all values of e, — oo < e < +00. More specifically, 
collimated solutions are found to the left of the thick solid 
branch for k > (mainly for e < 0) and to the right of 
the thick solid branch for k < 0. 

All branches converge for positive e and small Goo to- 
wards the thick solid line. Physically this corresponds to 
the limit (see Eqs. 4.4b and 3.17) where kV^/V.^ becomes 
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(b) and (d) of Figs. 1-3. 



negligible and the magnetocentrifugal coUimation arises 
because of the central EMR as in Paper III (e' « e). 
Thus we can compare Figs. la-3a to the related curve 
Knoo/2A2 = of Figs. 4 and 5b. As Hoc = 0, Goo, M^ 
and l^c>o/V^* depend only on e', namely 



M^ 



x/2 



2A2 



Gi 



4A2 



2\^/e' - 1 
4A2/e' - 1 



(4.5) 

This is identical to Eqs. (5.13) and (5.14) in Paper III with 
e' replacing e. 

However, even with e < 0, we have e' > provided 
that kV^/V^ is larger than |e|. In fact with Moo > Goo 
(^ Kx) > K), the ratio V^/V^ is large and positive and 
can compensate all negative values of e even for small val- 
ues of K. In this way, even as e — > —00, e' > and 
collimated solutions are obtained, albeit with rather large 
radii Goo- In physical terms, this corresponds to a sit- 
uation where the central source is an IMR and cannot 
collimate the flow through magnetic processes alone. Nev- 
ertheless, first it is possible that the conversion of thermal 
energy into kinetic energy is very efficient M^/C^ ^ 1. 
Second, it can be more efficient on a non polar streamline 
than on the polar one if there is more thermal energy in 
the nonpolar streamlines than in the polar one (k > 0). In 
this case then as the flow expands it will build up pressure 
gradients that will force the lines to bent towards the axis. 
Once the thermal energy is converted into kinetic energy 
and the pressure is becoming negligible the magnetocen- 
trifugal forces will dominate. However the coUimation is 
obviously less efficient (larger Goo) than that produced by 
a central EMR (Fig. la). Therefore, this excess of ther- 
mal energy induces the collimated character of the solu- 
tion through the energy integral, even though the corre- 
sponding pressure force fp does not enter directly into the 
asymptotic force balance condition. The term "magneto- 
centrifugal confinement" can be used for all {e/2X'^, k/2X'^) 
values of Figs. la-3a. 

For K < 0, in order to keep e' > 0, larger values of e 
are required, in comparison to the simple k = case. In 
physical terms this is so because now there is a deficit of 
thermal energy along the nonpolar streamlines in compar- 
ison to the polar one and therefore the star has to be a 
more efficient magnetic rotator (larger values of e > 0) in 
order to have a collimated outflow. This trend is shown 
by the grey branches to the right of the thick solid branch 
corresponding to k = (Figs. la-3a). 

Note that for k > 0, for each value of e/2X'^ there ex- 
ists a single value of Goo , Afoo and T4o • On the other hand 
for K < 0, more values of Goo, -Woo and T4o correspond 
to the same value of e. However the lower value of Goo is 
practically coincident with that of k = 0: in other words, 
thermal decoUimation is negligible, similarly to the case 
K = and e' « e. The upper values of the branches corre- 
spond to efficient conversion of thermal energy into polar 
acceleration but with a deficit along non polar lines. It in- 
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Fig. 5. Comparison of tlie plots of tlie asymptotic cylindrical radius normalized to the cylindrical radius at the Alfven surface 
Goo vs. e/2A^ (panel a) and vs. e'/2A^ (panel b). Each curves on the left panel (a) are drawn for a constant value of k,/2\^ 
which labels the curve. All these curves correspond to a unique curve on the right panel (b) because they all have the same 
value Knoo/2A^ = 0. We see that for each value of e /2\^ , we can find the same value for Goo for a pair of values of e/2A^ and 
«:/2A^ which have approximately the same magnitude but opposite sign (points 1 and 2). 



duces a more drastic conversion of the Poynting flux into 
acceleration along non polar lines thus reducing the effi- 
ciency of the collimation (Fig. 1, Goo larger) despite the 
fact that the central object is an EMR (e > 0). 

This is illustrated in Fig. 5 where the asymptotic cylin- 
drical radius Goo is plotted versus e/2A^ and e'/2A'^ for 
various k/2A^ and Hoc = 0. The higher value of Goo cor- 
responds to a higher conversion of thermal energy that 
decoUimates the wind which is balanced by a strong EMR 
as we explained. For some value of e' = e + nVi^/V^ , we 
can find the same value for Goo for a pair of values of e 
and K, which have approximately the same magnitude but 
opposite sign (see Fig. 5, points 1 and 2). 

Note that a situation with IIoo = corresponds to 
very specific initial conditions of integration that may not 
be easily fulfilled for a cylindrically coUimated flow. In 
particular IIoo = would imply Too = and the existence 
of some efficient cooling. 

4.3. Magnetocentrifugal jets with fc ~ |/b| 

In addition to the case wherein IIoo = shown in Figs. 
la-3a, where we have an exact magnetocentrifugal equilib- 
rium everywhere, approximate magnetocentrifugal equi- 
librium conditions also exist for IIoo t^ 0. This is shown in 
Figs. Ib-d, 2b-d, 3b-d and 4, in a region adjacent to the 
thick solid curve obtained for k = 0, between the dashed 
and the dotted lines. 



4.3.1. Underpressured jets, k > 

Following a branch of k > on the left side of the limiting 
curve K — 0, we see from the plots of the various forces 
shown in Fig. 6, that as e/2A^ decreases the magnetic con- 



finement of the jet is replaced by a pressure confinement, 
as expected. This transition from magnetic confinement to 
pressure confinement can be found by writing fp = fs, 
or equivalently, |/b| = /c/2 which gives 



G„ — 



M'i 



2Mo 



l + 2Mo 



(4.6) 



By inserting this relation in Eq. (4.4b), we obtain the dot- 
ted line of Figs. 1-4. Thus, for some finite value of IIoo, for 
each positive value of k/2A^, there exists a single value of 
e/2A^ = ep_s/2A^ located on the dotted line where there 
is an equal contribution by the magnetic and gas pressure 
forces in confining the jet against the outwards inertial 
(centrifugal) force. On the left side of the dotted line the 
jet enters the regime of gas pressure confinement, which 
we shall further discuss in the next subsection. 

We note that the limit between the two confinement 
regimes is close to the maximum of the cylindrical radius 
Goo as a function of e/2A^, for each value of k/2A^ (Figs. 
lb-Id). The limit is also very close to the minimum value 
that e'/2A^ can achieve for a given value of k/2A^ (Fig. 
4q, thick solid lines). On the other hand. Moo and Foo are 
monotonic functions of e/2A^. 

This maximum radius of the jet can be calculated for- 
mally from Eq. (4.4b), 



G 



Mi 



4M4 - 3M2 + 1 



(4.7) 

which can be combined with Eq. (4.4a) to give the radius 
as a function of the parameters e/2A2, k/2\^ and IIoo- 

It is interesting to note that, in the limit of large Alfven 
numbers (Afoo ^ 1) the two values of Goo given by Eqs. 
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(4.6) and (4.7) coincide and to first order we have 



G: 



(. 



V SkHo 



(4.8a) 



Then the asymptotic velocity along the polar axis is (see 
Eq. 3.4c) 






/ 8A2 



VkHo 



1/3 



A = 



fi(pole)r. 



(4.8b) 



The terminal speed has the same dependence on the di- 
mensionlcss rotational speed A with Michel's minimum en- 
ergy solution for cold magnetic rotators wherein V^oo/^* == 
\^/^ . In the present case however, the asymptotic speed is 
enhanced by the factor [8/Knoo]"'^^^. For small values of k 
and Hoc of the order of 1 , this is indeed a rather large en- 
hancement. This increased terminal speed simply reflects 
the transformation to asymptotic kinetic energy of the en- 
thalpy and added thermal energies. Note also that when 
Hoc — > 0, Ko — > oo. This is expected and the situa- 
tion is similar to the radial outflow studied in Tsinganos 
& Trussoni (1991) where the terminal speed is 



^«6(A2lni?' 
4.3.2. Overpressured jets, k < 



1/3 



(4.9) 



The regime of magnetocentrifugal equilibrium extends 
also to the right of the thick solid line n ~Q and up to the 
dashed line in Figs. 1 and 2. Beyond this line the jet en- 
ters in the regime of gas pressure support which wc shall 
discuss later. Following a branch of k < on the right 
side of the limiting curve k = 0, Fig. 7 illustrates how the 
centrifugal force decreases and is progressively dominated 
by the pressure gradient. The transition from centrifugal 
support to pressure support can be estimated by writing 
fp = fc, or equivalently, |/c| = /b/2 which gives 



GL = M„ 



(4.10) 



This relation can be combined with Eq. (4.4b) to give the 
dashed limiting line (Figs. 1-4) where e/2A^ = ep_c/2A^. 
For Goo ^ 1 this limit coincides also approximately with 
the maximum of the velocity on an assumed k/2A^ branch 
(Fig. 3), given by 



v; 



oo 



1 



2G„ 



[Gto - 4GL + 5 



VGL - 8G6^ 



26G^ - 40G^ 



21 



(4.11) 



(in this domain the curves of the jet radius Goo and Alfven 
number Moo are monotonic with e/2A^). We remark that 
the curve of maximum velocity (dashed line) also corre- 
sponds to the minimum value of e'/2A^ for a given value 
of k/2A^, in the limit of large Alfven number (see Fig. 47, 
thin solid lines). 



For Goo ^ 1 we have from Eqs. (4.10): 

2A2 ^l/^ 






Ma 



kn. 



(4.12) 



i.e. the same scaling law with |k|/2A^ and IIoo holds for 
the maximum velocity as for the asymptotic velocity at 
maximum radius for n > (Eq. 4.8b). 

4.3.3. ep-B < e < £p-c 

In the intermediate region bounded by the two curves 
ep_B and ep-c, the jet is magnetocentrifugal. Solutions 
very close to the thick solid line correspond to an ef- 
ficient collimation by the EMR with negligible thermal 
contributions, as we already discussed. For small values 
of Hoc (e.g., Figs. lb-3b) this area is surrounded by so- 
lutions with important thermal energy conversion but 
small asymptotic pressure gradients similar to the ex- 
tended branches of Figs. la-3a. However, for larger IIoo 
the jets can be in magnetocentrifugal regime only for a 
narrow range of values of e/2A^, around the line k — 
(see Figs. lc-3c and ld-3d: the region between the dotted 
and dashed lines shrinks by increasing Hoc). 

4- 4- Pressure confined jets fc ~ \.fp\ 

The more negative e/2A^ becomes, the weaker are the 
magnetic pinching forces (the less efficient is the magnetic 
rotator). Thus, for e < ep_B (left of the dotted line on 
Figs. lb-Id, 2b-2d, 3b-3d and Aa) the magnetic pinching 
force has dropped to very small values in comparison to 
the gas pressure force fp, such that now fp alone confines 
the jet against the inertial force (Fig. 6). The situation is 
similar to the case of the prescribed streamlines studied in 
TTS97 (q > 0, in the upper branch of curves in the right 
panel of Fig. 1 in this paper). 

The asymptotic radius of the jet and its Alfven number 
are sensitive to the nature of the asymptotic confinement. 
They strongly depend on the value of IIoo (e.g.. Figs. 1- 
2). At the same time, the terminal speed Voo/V^ is almost 
independent of the value of the terminal gas pressure IIoo 
(e.g.. Fig. 3) for strongly negative values of e/2A^. This 
can be understood from Eq. (4.4b) where we see that in 
the limit of negative e/2A^ the first term of the right hand 
side dominates and thus the square of the terminal speed 
is roughly given by the ratio |e|/K. 

For e ^ ep-B we are basically entering the hydrody- 
namic regime studied in Paper I. In the limit of e ^ —00 
and finite asymptotic pressure, the jet is strongly pressure 
confined such that Goo(e -^ —00) < 1 (cf. Figs. lb-Id, 2b- 
2d, 3b-3d), i.e., the solution becomes unphysical. As in Pa- 
per I we find that the most physically interesting hydrody- 
namic solution is obtained for vanishing terminal pressure 
with conical asymptotics wherein Goo(e -^ —00) -^ 00 
(e.g.. Fig. la - 3a). 
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Fig. 6. Asymptotic transversal forces vs. e/2A^ for a positive value of k, k/2A^ = 0.001: centrifugal force {fc, solid line), 
magnetic force (/s, dashed line) and pressure gradient (/p, dotted line). The asymptotic pressure is IIoo ~ (panel a), — 0.01 
(panel b), =0.1 (panel c) and = 1 (panel d). The bullet marks the value of the centrifugal force fc for which fp — /b- 



We may see the hydrodynamical limit as the most ex- 
treme one. Nevertheless, even with a non vanishing mag- 
netic field, we note that the more efficient is the pressure 
confinement of the jet (the more negative is e/2A^) the 
larger is the gap between the value of the jet radius Goo 
obtained for Hoc ~ and the one obtained for IIoo ^ 
(Figs, la and lb), for a given value of k/2A^. The differ- 
ence is even larger when k/2X^ takes small values. In such 
a case the magnetocentrifugal forces are very weak (Fig. 
6) and the equilibrium is very sensitive to small changes 
in the pressure gradient. 



4-. 5. Pressure Supported Jets with fp w |/b| 

This last case occurs when the centrifugal forces are neg- 
ligible i.e., the jet is confined by magnetic forces and is 
supported by the gas pressure gradient (for e > ep-c, 
right side of dashed line, but only in Figs. Ib-d, 2b-d and 
4a). 



Now the inertial force has dropped to very small val- 
ues in comparison to the gas pressure gradient force such 
that fp alone supports the jet against the magnetic pinch- 
ing force (see Fig. 7). In this domain, for a given value of 
k/2A^, the centrifugal force exactly vanishes for Moo = 
Goo (Eq. 4.3a), when T4o = V*. It simply states that a 
jet with zero asymptotic centrifugal force has no net ac- 
celeration between the Alfven surface and infinity. At this 
particular point e/2A^ — (e/2A^)o where the asymptotic 
centrifugal force is exactly zero and Voo — V^,, we obtain 
from Eq. (4.4b), 



V2A2Jo^ 2A2 



1 



(4.13) 



Since k/2X'^ is usually rather small, it follows that 



{e/2y 



1. There, we can say that the jet is exactly 



supported by the pressure gas gradient alone. If we start 
at (e/2A2)o and move toward larger values of e/2X'^, along 
the branch k/2X'^ = const < 0, Moo and Goo increase, 
rotation changes sign (Sec Eq. A.3f) and the asymptotic 
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a) n = 




.c)n =0.1 



6 / 2>/ 



K/2?i = -0.001 




b) n =0.01 



K/2?i = -0.001 




- :d) n„ = 1 



6 / 2A' 




K/2?i = -0.001 



E / 2\^ 



€ / 2A^ 



Fig. 7. The same as Fig. 6, for a negative value of k, k/2A^ = —0.001 but the magnetic force (/s) corresponds to the sohd 
line, the centrifugal force (/c) to the dashed one and the pressure gradient (/p) to the dotted one. The bullet marks the value 
of the magnetic force fs for which fp = fc- 



velocity is less than the velocity at the Alfven point (Fig. 
3) which means that the outflow is decelerated though it 
remains superAlfvenic (we do not consider in the present 
analysis "breeze" solutions that are always subAlfvenic) . 
By starting at (e/2A^)o and moving in the opposite di- 
rection toward smaller values of Moo and Goo along the 
branch k/2\^ — const < 0, the inertial force remains neg- 
ligible in comparison to the gas pressure gradient, up to 
the dashed line wherein we enter the magnctoccntrifugal 
domain. 

5. Oscillations in the jet's width 

The previous analysis giving the asymptotic equilibrium 
for confined jets can be pushed one step further to the 
first order terms in the bending of the lines (Paper III, 
Vlahakis & Tsinganos 1997, henceforth VT97). We also 
assume there that the jet becomes asymptotically cylin- 
drical. An expansion of G and M can be made then to get 
an idea of the fluctuations that exist far from the region 



of the initial acceleration of the wind, 
G'^GUl + ^), 



(5.1a) 
(5.1b) 



where ^, /i <C 1. Conversely to Paper III, we must also 
expand the pressure as 



n = noo(i+p). 



(5.2) 



Thus we obtain the harmonic oscillator equation for the 
perturbed jet radius (sec Appendix C for details): 



^+(27rr,/Ao,c.)'i?-0, 



(5.3) 



where Aqsc. is the wavelength of the oscillations. We can 
write the wavelength of the oscillations in the form of 
VT97, Eq. (28), namely 



27rr, 



2A2 



(1 - Mir 



m(2M 



l)Gt 



M. 



4 1 



M2^(l-M2^) 



(5.4) 
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Fig. 8. Plots of the wavelength of the oscillations Aosc. in units of r,/A vs. e/2A^ for the same values of IIoo and k/2A^ as 
in Figs. 1-3. On the left of the dotted line is the domain of pressure confined jets while on the right of the dashed line is the 
domain of pressure supported jets and in between is the domain of magnetocentrifugal jets. 



with in the present case 
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A2 



G^M^ + G^(1-2M^) 
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2(1 - M^ 
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(1-g: 
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(5.5) 



(1 - Ml,f 



Plots of the asymptotic -wavelength Aqsc. "vs. e/2A^ for four 
representative values of IIoo are sho-wn in Fig. 8. In each 
of these plots, the values of Aqsc. are plotted for the range 
of the pressure parameter k/2A^ as in Figs. 1 to 3. 

In the domain of magnetocentrifugal jets the -wave- 
length behavioiir is very similar to the one found in Paper 
III as expected.tl In the case of k = studied in Paper III, 



In Fig. 3 of Paper III, the plot of the wavelength in the 
region where Moo is close to one has to be corrected, due to 
the presence of a M^, at the denominator of Eq. (5.23) that 
should be replaced by a M^ . Fortunately this change does not 
affect the curve in the region of astrophysical interest, mainly 
for Moo > 1. 



Eq. (5.4) takes the very simple form (Vlahakis, private 
communication) 

4A2 



(27rr*/Ao 



Ml{M^^ - 1) 



(5.6) 



In the domain of pressure confined jets (-which corresponds 
al-ways to underpressured jets, k > 0), the wavelength of 
the oscillations behaves like the Alfven number despite of 
the decrease of the cylindrical radius. In particular, as the 
rotator slows down (e/2A^ decreasing) the wavelength in- 
creases very slowly. Nevertheless in the limit e/2A^ -^ — oo 
where we enter the hydrodynamic regime the wavelength 
must eventually diverge. The similarity of the curves of 
Figs. 2 and 8 on logarithmic scales, which reflects the sim- 
ilar behaviour of M^o and Aqsc , can be easily understood 
in the limit of large values of Moa and Kso/K • In this limit 
- which is by the way expected to be the case for most 
observed jets - we get that the wavelength is just propor- 
tional to the square of the Alfven number, Aqsc ~ M^ . 
Note also that, if the wavelength of the oscillations can be 
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related to the observed morphology of the jets, we may 
have here an indirect estimate of the magnitude of the 
poloidal Alfven number. 

For overpressured jets (k < 0), the behaviour of Aqsc. 
is similarly following the increase of the Alfven number 
as the pressure becomes more important, but as it enters 
the domain of pressure supported jets after the maximum 
velocity, the oscillations disappear. In fact in this last case 
we can see that the pressure gradient supporting the jet 
cannot restore equilibrium against the confining magnetic 
pinch. This could imply that the solutions of this class are 
unstable. 

We must notice finally that oscillations in coUimated 
winds are quite a general result, not restricted to our 
class of meridional self-similar solutions. In fact oscillating 
structures have been found not only in other self-similar 
flows (Chan & Henriksen 1980, Bacciotti & Chiuderi 1992, 
Del Zanna & Chiuderi 1996, Contopoulos & Lovelace 1994, 
Contopoulos 1995, VT98), but also in more general anal- 
yses of axisymmetric outflows (see e.g. Pelletier & Pudritz 
1992). 

6. Asymptotic equilibrium for non collimated flows 

6.1. Existence oj asymptotically non collimated flows 

If the flow is not cylindrical asymptotically, then Fao must 
take a value in the interval < F^o < 2, with F^o = 
corresponding to conical asymptotics. Assuming that this 
value of Foo and the corresponding shape may be achieved 
rather slowly, at large distances R from the central object 
the analytic expression of F{R) can be written as 



F{R ^ ^) = F« 



c 
h^R 



E 



R^' 



(6.1) 



where c and c„ are constants. As R — > oo the dominant 
term is the logarithmic one and thus we may keep only 
this term in the expansion, such that from Eq. (2.6) and 
(3.4c) we obtain 



G2(i? ^oo) = 



M'^iR 



R 



2-F„ 



./coln^(i?)' 



(6.2a) 



(6.2b) 



' fooln{R) K 

where /oo is a constant (TTS97, Paper III). By substitut- 
ing these expressions in the definition of the integral e of 
Eq. (3.11) and keeping the dominant terms we get 

1 y2 
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1 



T4o 



(6.3) 



Assuming k =/= 0, Eq. (6.3) shows clearly a rather gen- 
eral result: a diverging asymptotic velocity is inconsistent 



with the constancy ofe. The cylindrical radius and Alfven 
number may be unbounded, but the terminal flow speed 
cannot, no matter what is the exact value of F^c. This is 
also in agreement with the fact that the fraction of the 
heating term converted to kinetic energy and ex k cannot 
diverge (see Sec. 3.1). 

Further insight in the general behaviour of the asymp- 
totics can be gained by considering the dominant terms in 
the transverse momentum equation which gives a relation 
for the pressure (see Eq. A. 4c in Appendix A), 



Kn(R 



2A2 



= /o 



If hi%R)_ 



(6.4) 



The first conclusion from this relation is that the asymp- 
totic pressure must vanish because the r.h.s. term (due 
to the magnetocentrifugal terms) always vanishes for non 
cylindrically collimated outflows. 

A second conclusion is that, independently of the value 
of e, underpressured Bows (k > 0) must be cylindrically 
collimated, as found in Sec. 4. This can be explained tak- 
ing into account that, when the streamlines try to ex- 
pand, then pinching by both magnetic forces and pressure 
gradient will dominate over all other forces (Eq. 6.4). To 
maintain the forces equilibrium a strong bending of the 
lines towards the axis is required, RdF/dR > 0, such that 
the system must relax towards a collimated configuration. 
Then only overpressured outflows (k < 0) may be non 
collimated. 

6.2. Asymptotically paraboloidal or radial overpressured 
flows 

Assuming k < 0, we discuss separately the case of 
paraboloidal {Foo > 0) and radial asymptotics (Fqo — 0). 

6.2.1. Paraboloidal asymptotics, Foo > 

Independently on the value of c, Eq. (6.3) further simplifies 
as follows 

2A2 ^2A2 K2 Voo ' 



(6.5) 



which implies that paraboloidal shapes can exist only for 
e > 0. In other words, only overpressured outflows from 
an efficient magnetic rotator may eventually achieve a 
paraboloidal shape. This can be physically understood be- 
cause such a structure implies some coUimation that can 
be achieved only through magnetic forces in the present 
case. 

We remark that in the case k = e = Eq. (6.5) can be 
fulfilled only for diverging asymptotic velocity (Paper III) . 
However, as we have seen before, such a case would require 
an infinite heating rate, so that such kind of solutions 
should be considered as unphysical. 
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6.2.2. Radial asymptotics, Foo = 

li c =/= 0, Eq. (6.5) and the previous remarks still hold. For 
c = 0, Eq. (6.3) becomes 
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(6.6) 



and radial asymptotics can exist a priori for both negative 
and positive values of e. 

This simple analysis shows that overpressured flows 
{k < 0) can attain only radial streamlines if they are IMR 
(e < 0), while they can have cylindrical, paraboloidal or 
radial asymptotics if they are associated with an EMR 
(e>0). ^ 

We point out finally a common feature for winds with 
parabolic or radial asymptotics. From Eq. (6.4) we see that 
temperature goes to a constant value along each fieldline 



p 2 1 



5a 



-T)h '"' 



i.e. all uncoUiniated solutions are isothermal asymptoti- 
cally. This is consistent with the results of Tsinganos & 
Trussoni (1991), who analyzed solutions with prescribed 
radial streamlines [F{R) = 0]. More in general, we can 
expect that in non coUimated outflows some heating is al- 
ways necessary in the asymptotic regions. Conversely in 
a cylindrical jet, where the pressure and density are con- 
stant, the temperature is also constant but the heating 
rate in the flow vanishes (Eq. 3.1b,c). 

7. Discussion and astrophysical implications 

7.1. Summary of the main results 

We have presented here the asymptotic properties of su- 
perAlfvenic outflows which are self-similar in the merid- 
ional direction. The terminal Alfven number. Moo, the 
dimensionless asymptotic radius of the jet, C?oo and ve- 
locity, Foo/K depend only on three parameters (see Eqs. 
4.4). Besides the terminal pressure IIoo, the two other cru- 
cial parameters are: 

• k/2A^, connecting the radial and longitudinal compo- 
nents of the gradient of the gas pressure. Thus, the outflow 
can be either overpressured {k < 0), or underpressured 
{k > 0) with respect to the rotational axis. 

• e/2A^, which measures the magnetic contribution to 
the collimation of the outflow. Thus, we may divide the 
sources of outflows into two broad classes : Efficient Mag- 
netic Rotators (EMR) corresponding to positive values of 
e and a strong magnetic contribution to collimation, and 
Inefficient Magnetic Rotators (IMR) which have negative 
e and canjXoUimate outflows only with the help of the gas 
pressure. El 



■^ This does not imply that coUimated jets from IMR are 
asymptotically pressure confined. They may be magnetically 
confined at infinity if the terminal pressure IIoo is very small. 



The absolute criterion for cylindrical collimation is 
given by the sign of a combination of the two parameters 
k/2X^ and e/2X^. This new parameter 
• e' /2X^ is related to the variation across the streamlines 
of the various energy contributions which govern the flow 
dynamics (magnetocentrifugal, thermal, etc). Thus, a pos- 
itive value of e' provides cylindrical asymptotics while neg- 
ative values are required for having uncoUimated flows. 



Pressure confined 1 



IMR 



No Confinement 



Kin. 



Magnetically Confined 



EMR 



a I IX 



-- CYLINDRICAL Asymptotics 

= RADIAL (i.e. CONICAL) Asymptotics 



P = PARABOLOIDAL Asymptotics 



Fig. 9. A sketch of the various confinement regimes and 
asymptotical shapes of the fiow in the plane [k/2\^, e/2}?\. 
The pressure confined regime corresponds to the area filled 
with dark grey while the magnetically confined regime corre- 
sponds to the light grey filled area. In the remaining area only 
uncoUimated solutions are found. 



Despite this simple criterion, the asymptotic behaviour 
of the outflow still depends on the value of each parameter 
taken separately as shown in Fig. 9 [see also Figs. (1-3)]. 
• Underpressured outflows . For k > the wind always ob- 
tains cylindrical asymptotics. The jet is supported by the 
centrifugal force, while collimation can be provided either 
by the magnetic pinching or the gas pressure, depending 
on e and the value of the asymptotic gas pressure lloo. 
We repeat however that for IMR the state of asymptotic 



However in such jets the gas pressure always plays a crucial 
role in the achievement of the final collimation through strong 
pinching gradients in the intermediate region between the base 
and infinity. 
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magnetic confinement can be achieved only through the 
strong pressure gradients occurring between the base of 
the flow and infinity. If IIoo does not vanish, the jet radius 
has a maximum when moving from the magnetic to the 
thermal regime (by reducing e/2A^). 

• Overpressured outflows . When n < the jet can be con- 
fined by the magnetic pinch only, and is supported either 
by the centrifugal force or by the thermal pressure gradi- 
ent. If Hoc does not vanish, the jet terminal velocity has a 
maximum when moving from the centrifugal to the ther- 
mal regime. Moreover, in overpressured outflows, cylin- 
drical configurations are attained only for values of e/2A^ 
higher than some threshold depending on the pressure pa- 
rameters, IIoo and K. Below this positive threshold value 
the outflow reaches conical or paraboloidal asymptotics if 
e > (EMR) or only purely conical asymptotics if e < 
(IMR). 

Cylindrical coUimation seems to be quite a natural 
end state for superAlfvcnic outflows with non vanishing 
asymptotic pressure, as also found by other studies based 
on the radial self-similar approach (Li 1995, 1996, Ferreira 
1997, Ostriker 1997), or on the full asymptotic treatment 
of the MHD equations (Heyvaerts & Norman 1989). How- 
ever in the present case the coUimation can be provided 
not only by the magnetic pinch, but also by the thermal 
pressure gradient. This is consistent with our self-similar 
scenario, suitable to model winds close to their rotational 
axis, where the thermal effects are essential to drive the 
outflow. We also point out that our results are consistent 
with those of TTS97, where again different coUimation 
regimes can be found. Finally, we remark that cylindrically 
coUimated streamlines most of times undergo oscillations 
with different wavelengths. 

7.2. Astrophysical application 

As in Paper III and TTS97, the present results could 
be particularly suitable to model the physical properties 
of coUimated outflows associated with Young Stellar Ob- 
jects (YSO). However, since here we have analyzed only 
the asymptotic properties of winds, we will shortly dis- 
cuss only a simple possible scenario based on the physical 
difference between EMR and IMR. 

Let us consider a rapidly rotating magnetized proto- 
star at the beginning of its evolution. In such conditions 
this object could be considered as an EMR, with a well col- 
limated, magnetically confined jet. At the same time the 
pressure inhomogeneity \k\ and the asymptotic pressure 
may take rather high values due to the inhomogeneous 
and anisotropic environment in which the jet is found. In 
the early phases of stellar evolution the outflow extracts 
quite efficiently angular momentum from the protostar, 
reducing its spinning rate. From the point of view of our 
model this means that the system moves from the state of 
an EMR to that of an IMR as it lowers the value of e. Of 
course the details of the evolution may be more compli- 



cated due to the feeding of the wind by the surrounding 
accretion disk. Nevertheless the net end result should be 
a decrease of the spinning rate and subsequently of the 
efficiency of the magnetic rotator (Bouvier et al. 1997). 

If the jet is initially underpressured, for example being 
embedded in a dense molecular cloud, as it approaches 
the regime of the IMR a gradual widening of its radius 
is expected. This widening may reach a maximum value, 
wherein the magnetic and thermal contributions to the 
confinement are comparable. At this stage it is reasonable 
to expect a decrease of the asymptotic pressure because 
of this widening. At the same time, we may have a more 
homogeneous flow, i.e., effectively a decrease of n. Thus, 
in the subsequent evolution of the central young star, we 
expect that the jet continues increasing its radius. Accord- 
ing to this scenario, we can imagine an evolutionary track 
of the outflow along the maxima from Figs. Id to Fig. lb, 
wherein there is an equal contribution by the thermal and 
magnetic forces confining the plasma. In this sequence, 
the terminal velocity always increases. Figs. 3d to 3b, as 
K and the asymptotic pressure H^o get lower and lower 
values. Finally, the outfiow becomes a loosely coUimated 
wind. 

In the above scenario, it is essential to have a decrease 
of K and Hoo during the evolution. Because, otherwise, 
we see that an IMR (e < 0) together with a strongly un- 
derpressured plasma (high ulloo) have the result to over 
coUimate the wind with an asymptotic radius comparable 
to its Alfven radius. This result would be somehow in con- 
tradiction with the observed radii of outfiows from YSOs 
which are expected to be much larger than the Alfven ra- 
dius. Instead, we propose that both thermal and magnetic 
confinement decrease simultaneously during the evolution 
of the central source. 

However, jets from planetary nebulae (PN) may 
present a totally different situation where the primary 
source of confinement is a strong pressure gradient, 
kHoo 3> 0, associated with a source which is a very in- 
efficient magnetic rotator, although with a non zero mag- 
netic field. The terminal radii of jets from PN are indeed 
observed to be rather small after some huge initial widen- 
ing (Frank 1998). We also note that our analysis favours 
a hydrodynamical origin of jets from PN similarly to the 
GWBB model (cf. Mellema & Frank 1997, Frank 1998) 
and contrarily to a pure magnetic origin of the refocusing 
of the wind. 

The analysis of overpressured jets is more complex at 
first glance as different regimes are possible for the same e. 
Let us first consider that the jet is initially quite narrow, 
centrifugally supported and originates from an EMR; i.e., 
it is on the lower branches of the thin grey lines (k < 0) of 
Figs. 1-3. As in the previous case, the rotation slows down 
(e decreases), the outflow rapidly opens and the velocity 
rises. However, below a threshold value of e/2A^ the flow 
becomes uncollimated (see Fig. 9). 
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If, conversely, at the beginning the jet is pressure sup- 
ported (i.e., on the upper branches of the dotted hues in 
Figs. Ib-d, 2b-d), the behaviour is quite ambiguous: a re- 
duction or an increase of the jet radius and velocity crit- 
ically depends on the initial conditions {\k\, asymptotic 
pressure, jet radius) when the star starts slowing down 
its rotational rate. Then it is difficult to model the possi- 
ble evolution of a pressure supported jet. But, as we said 
previously, the absence of oscillations in this region may 
indicate that such equilibria are in fact unstable and never 
attained practically. 

It is evident from the above that the possible outflow 
evolution is critically related to its physical conditions, 
namely if it is underpressured or overpressurcd. In which 
of these regimes the outflow can be found depends on the 
detailed history of the wind. We remind that the thermo- 
dynamic conditions across the jet in its asymptotic regime 
depend also on the assumed structure of density (S) and 
the intensity of the gravitational field {ly) . These two pa- 
rameters do not enter in the present analysis, but are es- 
sential for the thermal acceleration of the wind (see Paper 
III and TTS97). Therefore the next step in the present 
analysis is to make a careful parametric study of the nu- 
merical solutions, solving the set of Eqs. (A. 4). This is 
also demanded in order to make a detailed comparison 
with observational data and implies that we construct so- 
lutions that connect the base of the flow with the super- 
Alfvenic region fulfilling the regularity conditions at the 
critical points (Sauty et al., in preparation). In particular, 
it will be crucial to see whether or not all the asymptotic 
regimes presented here can be attained. 



case the set of the MHD equations lead to a closed system 
whose treatment does not require any further 'a priori' 
assumption as the relation between the components of the 
pressure gradient (as in Paper III and in the present study) 
or the prescription of the streamline shape (as in TT97). 
Finally, we point out once more that, contrary to the 
radial self-similar studies, the role of the thermal struc- 
ture of the flow, which is related to the details of pro- 
cesses of input/output of heating in the gas, is essential in 
our model and not only in accelerating the flow but also 
in constructing its global shape. The problem of the ener- 
getic behaviour of astrophysical plasmas in different astro- 
physical contexts (solar and stellar coronae, YSO, AGN) 
is however still open. 
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A. Appendix A 

The classical equations of ideal MHD steady flows are 

V-B = V-(py) = Vx (y xB) = 0, (A.la) 

p{V ■V)V -{B-V)B/4Tr 

^-\/{p + By8TT)-g^er, (A.lb) 



7. 3. Future directions of study 

The present results have been obtained in the framework 
of a self-similar treatment of the axisymmetric MHD equa- 
tions. This implies some 'a priori' constraints on the struc- 
ture of the solutions, that we summarize in the following. 

The surfaces with the same Alfvenic number are spher- 
ical [M{R, 9) = M{R)], with the velocity vanishing on the 
equatorial plane. Furthermore the 9— and r— components 
of the gradient of the gas pressure are linearly related. 
These assumptions are not too constraining if we consider 
our model as suitable to describe the physical properties 
of the flow around the rotational axis. We remind also 
that our results are consistent with those found in TTS97, 
where the two components of VP were unrelated in a wind 
with prescribed cylindrical asymptotics. 

Such limits of the present treatment could anyway be 
overcome by a different scaling of the physical variables 
with the colatitude. It has been shown in VT98 that the 
assumptions of Eqs. (2.3) are just a particular case of a 
more general class of solutions, with no vanishing velocity 
on the equator and with a more complex expression for 
the pressure (another particular case is the one studied 
in Lima et al. 1996; see also Vlahakis 1998). In such a 



pV ■ [\/h 



VP\ 
P J 



(A.lc) 



where h is the enthalpy of the perfect gas, q is the local vol- 
umetric heating rate including true heating and cooling, Q 
is the gravitational constant, A4 is the mass of the central 
object and the other symbols have their usual meaning. 

Under the assumption of steady state and axisymme- 
try, the existence of free integrals as defined in the main 
text gives the usual forms for the poloidal (p) components 
of the magnetic field B and the velocity V, using spheri- 
cal (r, 9, if) or cylindrical {zu, ip, z) coordinates (for details 
see Tsinganos 1982) 



±jp X e^o , 



*4 

V = B 

" 47rp^' 



while the toroidal components are 

L{A)-^A 1 - n72 n{A)/L{A) 

7 1-M2 

L{A) m^n{A)/L{A)~M^ 



B^p — — - 



V^ = 



(A.2a) 



(A.2b) 



(A.2c) 



VJ 1 - M2 

where M is the poloidal Alfven Mach number as defined 
in Eq. (2.1). 
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Considering our assumptions, Eqs. (2.3), wc get that B. Appendix B 
the components of the velocity and magnetic field reduce 
in our model to 



B.l. On the variation oj the specific energy 



B,, 



Br — B* 

Be - -B, 
A 



1 



cos 6* . 



G2(i?) 
1 F{R) 



-B, 



G^{R) 2 
1-G^{R) 



sin0 , 



R sin 6 , 



Vr = V, 



Ve 



F^ = K 



X 



G2(i?)l-M2(i?) 
M'^{R) cos 61 

GHR) y^T+Sa{R,e) 
M^{R)F{R) sine 



GHR) 2 ^l + Sa{R,0) 
G^{R)-M^{R) i? sin 61 



(A.3a) 
(A.3b) 
(A.3c) 
(A.3d) 
(A.3e) 
(A.3f) 



G2(i?) l-M^R) ^l + SaiR,i 

These last equations, together with Eqs. (2.3), can be com- 
bined with the poloidal components of the momentum 
equation, Eq. (A. lb), to give three independent equations 
for four unknowns n(i?), F{R), G{R) and M'^{R). The 
system is closed with Eq. (2.6) for G. Two of these equa- 
tions arise from momentum-balance in the radial direc- 
tion, while the third one from momentum-balance in the 
meridional direction. Thus, we obtain 



dn 



2 
G* 



dM2 



dR 



-^{F-2 

i?2 ^ 



kFRH- kR^ — 
dR 



m^r^ 

2 dAf2 Sv^ 



= 0, (A.4a) 



F 



2i?G2 



R- 



,dF 
'dR 



F^ - F-2 



G2 dR M^G^ 
m2 r 



RG^ 
2\^R 1 - G2 



F2 

- + F-4 



2A2j? (M2 - G2)' 

'G2M2 (1-M2)2 G2~ (1 - M2)2 

2A2i?2 (1 _ 6*2)2 ^^2 



G2 (1-M2)3 dR 



{F-2G^) 
= 0, (A.4b) 



F dM2 



A2j?2 

(1-M2)2 



+ KR^G^n + i 


dF ^ 
^di?+^ 


M2 

2 


r dF f2 1 

YdR^ 2 -^\ 




r(M2- 

M 


^.2)2 ^ 
2 


^(1 - G^ 


f 



F^ - F - 2 



= . (A.4c) 



In Sec. 5 (osciUations of cylindrical jets), the first order 
expansion scheme of the previous momentum equations 
amounts to saying that we have the following expression 
for the force balance across the fieldlines 

p(Fp.v)yp - -^(v X Bp) X Bp 

47r 



B^p 

Attw 



V{zuBp 



— —vzu - 



VP. 



(A.5) 



In Eq. (3.3b), we find successively five terms which corre- 
spond to the variation, in units of the volumetric energy, 
of the magnetic rotator between any streamline (a) and 
the polar axis (pole) of 
(i) the poloidal kinetic energy, 



RHP 



El 

T 



1 



2A2 



p{R,a)L{a)n{a) 



1 



1 



X -piR,a)V;{R,a) - -p(ii',pole)y/(i?,pole) 

(ii) the gravitational energy, 

GM 



(B.l) 



5iy' _ 2A2 

"r" ^ ^ p{R,a)L{a)n{a)^' 



[p{R,a) -p(i?,pole)] , 

(B.2) 
(iii) the azimuthal kinetic energy (= along the polar 

axis) , 



A^ r 
G2 



M2-G^ 



n 2 



1- M2 



2A2 



p{R,a)L{a)n{a) 



(iv) the Poynting flux (— along the polar axis). 



piR,a)V^{R,a) 
(B.3) 



2A^ 



1-G2 



1 - M2 



2A2 



p{R, a)L{a)Q{a) 



n{a) 



*A(a 



■zuBp{R, a) 



2A2 



p{R,a)L{a)Q{a) 
(v) the thermal content, 

r 



L{a)n{a) - m^{R, a)n^{a)] 



(B.4) 



r- 1 



-BM^- 



Q{R)dR 



2A2 



p{R^ oi)L{a)Q.{a) 



X < p{R, a)[h{R, a) 



el(«)] 



-p(i?, pole) [/i(i?, pole) - e^Jpole)] 



(B.5) 



In order to write Eqs. (3.7), we simply calculate the pre- 
vious terms (B.1-B.5) at the base of the flow ro 



-2A2 



5u^ 
Ro 


A2 

+ K 


\M^-Gl] 
[ 1 - M2 J 

r 

n If 


2 


r-i"°^'^ 


o 



(B.6) 



There the poloidal kinetic energy is negligible and conse- 
quently the term in Eq. (B.l) is zero. Moreover the ther- 
mal content reduces to the enthalpy such that in Eq. (B.5) 
only the enthalpic terms remain. Combining Eqs. (B.2) to 
(B.6), we get Eqs. (3.7). 



Sauty et al. Nonradial and nonpolytropic astrophysical outflows. IV. 



23 



B.2. On the Energetic definition of EMR and IMR 
Now using Eqs. (3.2) and (3.10) we may write 



E{a) 



i^^g^Ci + '*'^) + '^^ 



:VJ 



l + da 



(B.7) 



where we see, as stated in the main text, that e is the 
transverse variation of the volumetric energy once we have 
removed the thermal terms that linearly scale with factor 
K. Moreover note that from Eq. (3.4a) we get 



S _ /Oo(polc)[/i,o(polc) 
2A2 ~ Po(a)E-i 



E, 



G.oJ 



MR 



(B. 



which is the pending expression to Eq. (3.7a). Noting that 
A{ph) AP 



p(pole)/i(pole) P(pole) 
we combine Eqs. (3.7a), (B.8) and (3.10) to get 

A [po(£^Poynt.,o + Ey{,.o)] 



(B.9) 



e 
2A2 



Po{a)EMR 



Eg,o /Oq (pole) 
Emr Po{a) 



Apo 



APo 



Po(polc) Po (pole) 

The first term of the r.h.s. of this last equation simplifies 
to the two first terms in the numerator of Eq. (3.12a) as 
the Poynting flux and the rotational energy vanish along 
the pole. The second term of the r.h.s. of Eq. (B.IO) can 
be rewritten to give AEq in Eq. (3.12a). In the form pre- 
sented in Eq. (B.IO) it appears how the relative increase 
of the weight of the plasma can be partially compensated 
by the increasing of the thermal pressure gradient. In this 
form there is no contradiction with the use of the sym- 
bol A. Conversely, the equivalent expression used in Eq. 
(3.12b) may appear confusing if one does not remember 
that this is in fact the variation across the lines of the grav- 
itational energy that is not compensated by some thermal 
driving. Nevertheless, we prefer this last form for its com- 
pactness and because it emphasizes the role of the tem- 
perature. 

C. Appendix C 

From Eqs. (5.1) and (5.2) we may expand E to first order 
in Eq. (2.6), 

E = 2-m, (C.la) 

while the derivative of F can be also expanded at large R 
as 



dF 
OR 



= -t?-i?^ w -R^t}. 



(C.lb) 



From these equations we can expand the momentum equa- 
tions given in Appendix A (Eqs. A. 4) replacing the second 
one by the definition of e. We still assume in this Section 
that the fiow is asymptotically cylindrically coUimatcd. 
Thus we already have calculated the zeroth order equilib- 
rium in Sec. 4. We know that the asymptotic quantities in 



the flow arc uniquely determined by the values of e/2X'^, 
k/2X'^ and Boo- The first order terms in the transverse 
momentum equation (see also Eq. A. 5) give a relation be- 
tween ^, "d and /i 



?9 = 2kII^ 
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^-Ml 



-2\' 



2Gl 



-d- 



l-Gl 



2Ml 



1 - Ml. 



Gl 



{I-MIY \l-M^^ Ml 



'(1-MIYMI^_ 



(C.2) 



This can be combined with Eq. (3.11) that we also 
expand to first order - where again the zeroth order is 
Eq. (4.4b) ~ to get a second relation between p and i?. 
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0. 



(C.3) 



2(1 - Mir 

which is identical to Eq. (5.5). Expanding Eq. (A. 4a) we 
have a relation between the derivatives of 9, p and p, 



Mi 



p 



(-foo^^oo 



(A-^)-o, 



(C.4a) 



that can be integrated to give p as a function of 9 and p, 
assuming a vanishing constant of integration 



Mi 



(-foofl^oo 



(p-d)^O. 



(C.4b) 



By eliminating p and p in Eq. (C.2) using Eqs. (C.3)- 
(C.4b), we obtain Eqs. (5.3) and (5.4) in Sec. 5. 
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